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Statement of Significance 

QUESTION: Most land surface models (LSMs) currently used with general 
circulation models (GCMs) do not include adequate treatments of the lateral 
redistribution of subsurface moisture. This severely limits their production 
of realistic surface fluxes. How can this fundamental deficiency be addressed? 

APPROACH: Our paper summarizes a new strategy for modeling land sur- 
face processes. The hydrological catchment is defined as the basic land sur- 
face element, and within each catchment, the spatial heterogeneity of soil 
moisture is diagnosed from bulk moisture variables and the catchment’s to- 
pography. This spatial heterogeneity is an important control of evaporation 
and surface runoff. The paper alse describes the numerical framework to test 
the validity of this approach. 

SIGNIFICANCE AND IMPLICATIONS OF FINDINGS: The explicit con- 
sideration of subgrid-scale soil moisture heterogeneity and its effects on the 
surface energy and water budgets should lead to an improved simulation of 
land surface processes. This should in turn lead to more reliable simulations 
of climate and its sensitivity. 


RELATIONSHIP TO MTPE SCIENCE PLAN 

The improved simulation of land surface processes should benefit numerical 
climate studies. The strategy described in the paper thus has direct relevance 
to MTPE emphases on (1) seasonal-to-interannual climate prediction, (2) 
changes in long term climate, and (3) landcover and land use change. 
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Abstract. A new GCM-scale land surface modeling strat- 
egy that explicitly accounts for subgrid soil moisture vari- 
ability and its effects on evaporation and runoff is now be- 
ing explored. In a break from traditional modeling strate- 
gies, the continental surface is disaggregated into a mosaic 
of hydrological catchments, with boundaries that are not dic- 
tated by a regular grid but by topography. Within each catch- 
ment, the variability of soil moisture is deduced from TOP- 
MODEL equations with a special treatment of the unsatu- 
rated zone. This paper gives an overview of this new ap- 
proach and presents the general framework for its off-line 
evaluation over North-America. 


1 Introduction 

The realistic representation of land surface processes is criti- 
cal for the realistic simulation of the global hydrologic cycle 
and climate, as indicated by numerous sensitivity studies us- 
ing general circulation models (GCMs) (e.g. Garratt, 1993; 
Mintz, 1984). Land surface models (LSMs) have therefore 
increased in sophistication and realism over the last decade. 
An important effort has focused on the modeling of evap- 
oration, which controls both the surface energy and water 
budgets. In doing so, the vertical dimension of evaporation 
processes was largely emphasized, with detailed treatments, 
for instance of canopy structure and environmental stresses 
on surface conductances. 

The Mosaic model (Koster and Suarez, 1992) is represen- 
tative of the type of LSM that emphasizes evaporation pro- 
cesses. It includes a canopy interception reservoir and a 3- 
layer soil reservoir. Interception loss, bare soil evaporation 
and transpiration occur in parallel and are controlled by resis- 
tances that increase with environmental stresses, according 
to SiB formulations (Sellers et al., 1986). A distinct feature 
of the Mosaic model is its treatment of subgrid scale hetero- 
geneity in surface characteristics, which motivated its name. 
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The grid cells are subdivided into homogeneous sub-regions 
(the “tiles” of the mosaic), each containing a single vegeta- 
tion or bare soil type. A separate energy balance is calculated 
for each tile, and each tile maintains its own prognostic soil 
moisture and temperatures. But the Mosaic model, like most 
of the state-of-the-art LSMs, features overly crude parame- 
terizations of runoff and baseflow processes relative to the so- 
phistication of evapotranspiration, albedo, and other energy 
balance formulations (Koster and Milly, 1997). In particu- 
lar, the critical effects of horizontal soil moisture variability 
are rarely accounted for in current LSMs, and soil moisture 
is assumed to be uniform at the GCM scale. This is over- 
whelmingly unrealistic given the actual spatial scales of soil 
moisture variability, which are largely imposed by topogra- 
phy and precipitation heterogeneity. 

Most of the few attempts to describe soil moisture vari- 
ability at the GCM scale made use of statistical distribu- 
tions, either of precipitation rate or of soil properties (e.g. 
Duchame et al., 1998; Entekhabi and Eagleson, 1989). Fol- 
lowing recent advances in catchment models for the large 
scale (Famiglietti and Wood, 1994; Stieglitz et al., 1997), 
we are examining a new approach to the modeling of land- 
surface processes in GCMs. This approach recognizes the 
hydrological catchment as the fundamental land-surface unit, 
and calls for the disaggregation of the land-surface into a 
mosaic of hydrological catchments, with boundaries that are 
not dictated by a regular atmospheric grid but by topogra- 
phy. The second step of the new approach consists in the use 
within each catchment of a model that relates soil moisture 
variability to topography. A first incarnation of this catch- 
ment model is briefly outlined in section 2 and the construc- 
tion of the catchment coverage is outlined in section 3. Sec- 
tion 4 describes the framework for the validation of the new 
approach, which is discussed in section 5. 

2 Description of the catchment model 

It is widely recognized that a major control on soil moisture 
heterogeneity is topography (e.g. Beven and Kirkby, 1979; 
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Western and Grayson, 1998; Yeh andEltahir, 1998). Regions 
of local concavity tend to be zones of convergent flow, with a 
shallow water table and a high soil moisture content. Those 
regions are therefore characterized by high runoff (overland 
and baseflow) and by evaporation rates close to potential. In 
comparison, upland soils tend to be progressively drier, with 
a small contribution to saturation excess runoff (or Dunne 
runoff), and evapotranspiration rates significantly less than 
potential because of water stress. 

Beven and Kirkby (1979) developed the simple physically- 
based TOPMODEL to describe the effect of topography on 
runoff generation and catchment hydrology. This model has 
two fundamental assumptions. The first one is that the wa- 
ter table shape is related to topography, and the second one 
states that at any point in the catchment, the water table depth 
results from quasi-steady state conditions. TOPMODEL can 
then produce the spatial distribution of the water table depth 
across the catchment, given the mean water table depth, the 
topography of the catchment and the vertical distribution of 
soil hydraulic conductivity. The latter is usually assumed to 
decrease exponentially with depth, following Beven (1984), 
but recent developments on TOPMODEL allow to use differ- 
ent transmissivity profiles (Ambroise et ah, 1996; Duan and 
Miller, 1997). 

Although TOPMODEL was first proposed for small to me- 
dium sized catchments in humid areas, it has been since adap- 
ted to a broader range of scales and hydro-climatic condi- 
tions (e.g. Sivapalan et al., 1987). The catchment model ex- 
amined in this paper is based on some ideas from Famiglietti 
and Wood (1991, 1994) to extend TOPMODEL concepts to 
the macroscale. 

Our modeling strategy will be described in detail in an up- 
coming publication. Briefly, we characterize the moisture 
state in the catchment with two bulk moisture variables. The 
first variable is the “catchment deficit”, which is defined as 
the average amount of water, per m 2 , that must be added to 
the catchment’s soil to bring the entire catchment to satura- 
tion. The calculation of this deficit is based on the distri- 
bution of water table depths derived from the TOPMODEL 
equations and on an assumption of equilibrium soil moisture 
profiles in the unsaturated zone. The second bulk variable, 
termed the “root zone excess”, represents the average amount 
of water, per m 2 , by which the root zone moisture is differ- 
ent from the value implied by the equilibrium profiles. This 
second variable is important for two reasons: (1) it allows 
the root zone moisture to respond quickly to storm events 
and evaporation, and (2) it essentially makes the catchment 
model equivalent to a more traditional LSM with vertical soil 
layers in regions of little topography, where TOPMODEL is 
known to be invalid. Transfers between the root zone excess 
and the equilibrium state proceed according to timescales 
that vary with the magnitudes of the variables. 

The bulk moisture variables are combined with character- 
istics of the topography to derive a distribution of root zone 
soil moisture, which in turn is used to separate the catch- 
ment into hydrologic regimes. This is illustrated in Figure 1. 
As shown in the figure, the transpiration, bare soil evapora- 


pdf = f(M,M x , topography) 
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Fig. 1 . Fractioning of the catchment into four hydrological regimes. WP 
designates the wilting point, STR a vegetation-dependent threshold between 
stressed and unstressed transpiration; SAT designates saturation, BSE bare 
soil evaporation, M the catchment deficit and M x the root zone excess and 
pdf stands for probability densitity function. 


tion, infiltration, and baseflow calculations are different for 
the different regimes. We thus treat explicitly the effects of 
subgrid soil moisture variability on the surface energy and 
water balances. 


3 Construction of the catchement coverage 

The US Geological Survey has developed a global digital el- 
evation model (DEM) at 30-arc-second resolution (approx- 
imately 1 kilometer). This DEM was treated over North- 
America with a geographical information system to compute 
local slopes, and drainage directions were then extracted ac- 
cording to a “steepest descent” algorithm (Jenson and Do- 
mingue, 1988). This allowed the construction of the North- 
American drainage network. Catchments have been delin- 
eated according to ordering rules based on network topology 
and drainage areas (Verdin and Jenson, 1996). The ordering 
rules are associated with a catchment coding system (Pfaf- 
stetter, 1989). One of the appeals of this system arises from 
the simple way to assess upstream-downstream relationships 
between the catchments from their code alone. 

This procedure resulted in a coverage of 5020 catchments 
over North- America, with an average catchment size of 3,640 
km 2 . Figure 2 displays a fraction of that coverage in the 
South-western part of North America, with overlain l°x 1° and 
4° x 5° grids, the latter being representative of GCM grids. 
The figure shows the mismatch between the catchment cov- 
erage and the grids, as well as the difference in resolution. 
Most of the l°x 1° grid cells hold several catchments, and 
the 4° x 5° grid hold many of them. Using the catchments 
as the fundamental hydrological units therefore allows us to 
account for a first level of GCM subgrid scale variability. 

Topography in TOPOMODEL, and in the new model, is 
described in terms of a topographic index, the distribution 
of which can be idealized by a three-parameter gamma dis- 
tribution (Sivapalan et al., 1987). The three parameters are 
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Fig. 2. Comparison of the catchment coverage with l°x 1° and 
4°x 5° grids. 


a function of the first three moments (mean, variance and 
skew), minimum and maximum of the actual distribution of 
topographic indices in the catchment. 

The computed slopes and drainage areas have been used 
to compute the topographic index at every GTOPO30 grid- 
point in North-America. These values are used to determine 
the distribution of the topographic index in every catchment 
and the statistics of those distributions (mean, variance, skew, 
minimum and maximum) are currently being converted for 
each catchment into topography-related parameters used by 
the catchment model. 

4 Validation strategy 

4.1 Off-line Method 

The accuracy of a LSM can only be determined through a 
detailed comparison of its products with observations. How- 
ever, the validation of a LSM coupled to an atmospheric 
model can be very difficult if the forcing by the atmospheric 
model is in error. For instance, simulated precipitation rates 
are notoriously inaccurate, and these errors may dominate 
the behavior of the LSM. This explains the widespread use 
of “off-line” simulations, for which the atmospheric forc- 
ing is derived from observations, in validation studies (e.g. 
Dirmeyer and Dolman, 1998;Henderson-Sellerset al., 1996; 
Shao et al., 1994). 

The ISLSCP Initiative I data set for 1987-88 (Sellers et ah, 
1996) contains all of the atmospheric and boundary condi- 
tions data sets needed to drive a LSM off-line at the 1° x 1° 
resolution across the globe. Global fields of precipitation, 
incoming long-wave and short-wave radiation at the surface, 
and near-surface temperature, humidity, pressure and wind 
speed, have various observational sources. They are pro- 
cessed with a data assimilation system in order to extrapo- 
late them spatially across the globe and to interpolate them 
to the 6-hour timescale. The ISLSCP Initiative I data set 


also provides with data sets of vegetation characteristics at 
the 1° x 1° resolution, as derived from satellite observations. 
The data include vegetation type as well as monthly vary- 
ing fields of leaf area index, roughness length and snow-free 
albedo. 

The general framework for the off-line integration of the 
new model with the ISLSCP Initiative I data set for 1987-88 
is based on the recommendations for the Global Soil Wet- 
ness Project (Dirmeyer and Dolman, 1998). In particular, 
to avoid non-equilibrated spin-up signal, the forcing corre- 
sponding to 1987 is repeated until every modeled catchment 
has converged to an equilibrium state. All the needed forc- 
ing fields have been interpolated from the 1° x 1° grid to the 
catchment space. The difficulty was to determine the inter- 
section between the catchments (defined as polygons) and the 
grid cells, in order to compute the forcing in one catchment as 
the weighted average of the forcing from all overlying grid- 
cells. This has been achieved with a triangulation algorithm, 
which can also be used to transform data from the catchment 
space to the grid space. Those transformations in the two di- 
rections are of critical importance given our eventual goal of 
coupling of the new model with a GCM. 

4.2 Validation of annual runoff rates 

The off-line framework limits the number of available vari- 
ables for validation, since many of commonly measured vari- 
ables are prescribed. Three potential variables for valida- 
tion are runoff, evaporation, and soil moisture. Runoff is 
a choice variable for validation purposes because it can be 
easily estimated from river discharge measurements, which 
have been collected for many decades in many areas. Also, 
river discharge naturally integrates runoff across large dis- 
tances, which makes its measurement more valuable at the 
large scale than point measurements of evaporation or soil 
moisture. 

4.2. 1 Runoff observations 

The original observations consist of annual discharge at a 
collection of stations scattered across the globe. These data 
were processed in the manner described by Koster et al. (1998). 
The suitability for validation purposes of a river basin (de- 
fined as the direct contributing area for one station if two or 
more stations lie on the same river network) is assessed ac- 
cording to two criteria: (1) discharge was accurately mea- 
sured during 1987-88, and (2) at least 30 rain-gauges per 
10 6 km 2 were present is the basin at that period. The lat- 
ter criterion is necessary because the accuracy of simulated 
runoff is conditioned by the accuracy of the prescribed pre- 
cipitation. Those criteria lead to the selection of 14 river 
basins in North-America. Figure 3 shows their location de- 
rived from the catchment coverage (as explained in the next 
paragraph). For instance, the biggest basin is part of the Mis- 
sissipi river basin, and three sub-basins are found along both 
Saint-Lawrence River. 
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Fig. 3. Location of the selected 1 4 river basins as derived from the catchment 
coverage. 

4.2.2 Aggregation of simulated runoff 

The catchment ordering scheme mentioned in section 3 de- 
scribes the upstream-downstream relationships between the 
catchments. It therefore allows us to identify the catchments 
that contribute to the discharge measured at the 14 selected 
stations. As mentioned before, Figure 3 shows the resulting 
location of the 14 selected river basins. The area of each 
unit catchment being known, it is straightforward to com- 
pute the area of each basin, as well as the weighted average 
of runoff within those basins. Figure 4 plots these area esti- 
mates against the corresponding areas from observation. The 
strong correlation between observed and estimated areas is 
an evidence of the accuracy of the drainage network extrac- 
tion from the DEM. 

4.2.3 An objective criteria for validation 

An important limitation of off-line validation is the lack of 
feed-back between land-surface and the atmospheric forcing. 
This is the price we pay to eliminate atmospheric modeling 
errors when validating LSMs, and it can result in an over- 
whelming influence of the prescribed atmospheric forcing 
on the simulated land surface fluxes. For instance, even the 
poorest FSM can produce realistically low evaporation rates 
if forced with realistically low precipitation rates. A critical 
problem in off-line validation is thus to quantify the relative 
control of prescribed atmospheric forcing on the simulated 
energy and water balances, to assess which part of the LSM 
response is really due to its own structure. 

In order to validate the new LSM in such a objective way, 
we will apply the simple criteria proposed by Koster et al. 

( 1 998). It is based on an equation proposed by Budyko (1974) 
that relates annual evaporation to annual precipitation and net 
radiation, without any consideration of the surface physics. 
In an off-line framework, annual runoff can be approximated 
by the difference between annual precipitation and annual 
evaporation, under the assumption that the inter-annual vari- 



Fig. 4. Comparison of the area of the 14 selected river basins from observa- 
tion and from the DEM. The solid line comes from linear regression, and r 
is the correlation coefficient. 


ations of soil moisture are negligible in comparison with the 
fluxes themselves. The equation then provides with an an- 
nual runoff estimate that can be seen as a “yardstick of suc- 
cess”: if the formulations of a LSM contribute to the realism 
of the fluxes, then the LSM’s runoff estimates should be bet- 
ter than the estimate from Budyko’s equation. This estimate 
has been computed for the 14 selected basins, using annual 
precipitation and radiation data interpolated in the catchment 
space (section 4.1). The runoff simulated with the new LSM 
will have to be significantly closer to the observed runoff than 
to that estimate, for the new approach to be validated at the 
annual timescale. 

The above validation strategy favors the annual timescale, 
but seasonal and diurnal cycles have to be realistically sim- 
ulated by LSMs for realism in the GCMs short-term dynam- 
ics. In many ways, land-surface physics is more important at 
short timescales (e.g. Chen et al., 1997; Koster et al., 1998). 
The runoff data presented in section 4.2.1 are available as 
monthly means, and we plan then to use those data to eval- 
uate the behavior of the new LSM at the seasonal timescale. 
Care will be taken to minimize problems related to the time 
lag induced by river routing and regulation. 


5 Summary and discussion 

Many studies have shown the importance of small scale vari- 
ability of hydrological processes in the average interactions 
between the land surface and the atmosphere. This impor- 
tance limits the ability of point-process models to represent 
large scale hydrological processes realistically. These con- 
siderations motivated the development of our new catchment- 
based land surface model for use in GCMs. 

This new model calls for the disaggregation of the land 
surface into a mosaic of hydrological catchment, determined 
through the analysis of a high resolution DEM. Within each 
catchment, a new LSM based on TOPMODEL accounts for 
topographically-driven soil moisture variability and its ef- 
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fects on evaporation and runoff. The model introduces the 
root zone excess prognostic variable to ensure that the lim- 
itations of TOPMODEL in regions of little to moderate to- 
pography are minimized. 

Validation is a critical element of model development, and 
the second part of the paper was devoted to the description 
of the off-line validation framework for our new modeling 
approach. The very preliminary results using that framework 
do not yet show a significant improvement in annual runoff 
predictions when going from the Mosaic model to the new 
LSM. One must however keep in mind that this new model 
is still under development; not all of the model’s intended 
features have yet been included. We plan for instance to de- 
vise a sensible scheme for allowing surface runoff to occur 
from the unsaturated areas of the catchment, and we will in- 
clude some representations of bedrock to limit the attainable 
depth of the water table. Both of these features should in 
particular increase runoff, which is currently biased low. The 
main point we want to stress is that the new catchment strat- 
egy avoids the inherent problems associated with the one- 
dimensional representation of subsurface soil moisture. Be- 
cause of this, we expect the model, once fully developed, to 
produce a more realistic simulation of land surface processes 
than is achievable with more traditional approaches. 
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